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^ ABSTRACT 

Context. To address the physics behind the soft X-ray light curve asymmetries in Cygnus X-3, a well-known microquasar. 

^ Aims. Observable effects of the jet close to the line-of-sight were investigated and interpreted within the frame of light curve physics. 

D Methods. The path of a hypothetical imprint of the jet, advected by the WR-wind, was computed and its crossing with the line- 

of-sight during the binary orbit determined. We explore the possibility that physically this 'imprint' is a formation of dense clumps 
triggered by jet bow shocks in the wind ("clumpy trail"). Models for X-ray continuum and emission line light curves were constructed 
using two absorbers: mass columns along the line-of-sight of i) the WR wind and ii) the clumpy trail, as seen from the compact star. 
These model light curves were compared with the observed ones from the RXTE/ ASM (continuum) and ChandrafHETG (emission 
lines). 

L ' Results. We show that the shapes of the Cygnus X-3 light curves can be explained by the two absorbers using the inclination and true 

M— j anomaly angles of the jet as derived in Dubus et al. (2010) from gamma-ray Fermi/LAT observations. The clumpy trail absorber is 

I I I much larger for the lines than for the continuum. We suggest that the clumpy trail is a mixture of equilibrium and hot (shock heated) 

• clumps. 

Conclusions. A possible way for studying jets in binary stars when the jet axis and the line-of-sight are close to each other is 
demonstrated. The X-ray continuum and emission line light curves of Cygnus X-3 can be explained by two absorbers: the WR 
Q companion wind plus an absorber lying in the jet path (clumpy trail). We propose that the clumpy trail absorber is due to dense 

^ clumps triggered by jet bow shocks. 

H— > 
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^ 1. Introduction 

> 

Cn Cygnus X-3 (4U 2030H-40, V1521 Cyg) is a high-mass X-ray 
On binary (HMXB) located at a distance of 9 kpc with a close bi- 
nary orbit (P = 4.8 hour; Hanson et al. 2000; Liu et al. |2007] |. 
The compact star is either a neutron star or a black hole and the 
p\| companion is most probably a WN5-7 type Wolf-Rayet (WR) 
T— I star (helium star; van Kerkwijk et al. 119921 11996i However, it 
is possible that the WR-phenomenon comes from the accretion 
disc wind, but to our knowledge no such detailed modeling for 
this source exists. Cygnus X-3 is also a source of relativistic jets 
• ^ (e.g. Mioduszewski l20(JTt Marti ISOOTb . thus associating Cygnus 
X-3 with the class of microquasars. 
H Massive winds are generally observed in Wolf-Rayet stars 
(Langer 1989; Crowther I2007I I. This wind produces the well- 
observed orbital modulation of radiation from the compact star 
by asymmetric absorption during the orbit along the line-of- 
sight. 

The system incUnation (/=the angle between the line-of-sight 
and the orbital axis) is unknown but probably small. Depending 
on inclination, the compact star mass is either large (3-10 M0 if 
i = 30 deg) or smaller (1 - 3 M© if / = 60 deg) (Vilhu et al. |2DU9l 
see their Fig. 10). By modeling the X-ray spectra in all states, 
Hjalmarsdotter et al. (2009) concluded that the compact star is a 
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massive black hole (30 M0) indicating a very small inclination. 
Koljonen et al. (2010) found six distinct X-ray states reminiscent 
of those seen in black hole transients, also supporting the black 
hole case. 

Hot star winds are clumpy rather than homogeneous, con- 
sisting of dense condensations inside a rare gas (see papers in 
Hamann et al. 2008 ). In Cygnus X-3 the clumps are photoionized 
by strong X-ray/UV radiation and are consequently highly radia- 
tive. Jets destroy clumps when colliding with them (Perucho & 
Bosch-Ramon 1201 11 1. However, the opposite may also be true, 
such as in star formation regions where clumps instead form 
when a shock front collides with interstellar clouds (van Breugel 
etal. 12004b . 

In this paper, we explore the possibility that a bow shock 
from the jet indeed triggers the formation of clumps and/or 
makes the existing clumps denser and that these are advected 
by the WR-wind. The result is a "clumpy trail", a volume filled 
by the clumps and advected by the wind. In reality the formation 
of such a trail is far from being demonstrated but could in prin- 
ciple be investigated by radiative hydrodynamics (see Perucho 
and Bosch-Ramon I20()8t . The continuum and line absorptions 
are sensitive to density, and in different ways in a photoionised 
medium. Hence, one can expect observable effects in X-rays if 
the line-of-sight from the compact star (X-ray source) crosses 
the clumpy trail region. 
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We use the parameters derived by Dubus et al. (2009) for the 
jet inclination and true anomaly angles, and compute the inter- 
section of the clumpy trail with the line-of-sight. Two different 
bow shock geometries are applied, and the continuum and emis- 
sion line light curves modeled and compared with the observed 
ones. 

2. Geometry 

Dubus et al. ( 2010b modeled the FermifLP^ gamma-ray orbital 
modulation of Cygnus X-3 by inverse Compton scattering of 
WR photons on the jet using two models for the compact star 
given in Szostek & Zdziarski (2008): a neutron star with high 
orbital inclination (60 deg) and a black hole with small inclina- 
tion (30 deg). Dubus et al. derived values for the jet inclination 
and true anomaly angles, and showed that the jet axis is close to 
the line-of-sight. We used the black hole alternative, supported 
by the spectral studies of Hjalmarsdotter et al. (2009), Vilhu et 
al. (2009), and Koljonen et al. (2010). The WR-wind is assumed 
to be spherically symmetric with density depending on the dis- 
tance r from the WR-companion center as r^~^ with y -2. 

The rectangular xyz -coordinate system was used where the 
WR star is in the center and the z-coordinate is parallel to the 
orbital axis ( see Fig. [T). The vector formalism and definitions 
given in Dubus et al. (2010) were used. The unit vectors of the 
jet axis and line-of-sight (observer) in the xy z-coordinate system 
are: 

Sje, - (cos0ysin0j, s,m@ jS,m(pj, cos(pj) 
Sobs = (0, -sin/, cos/). 

Here 0^ is the jet azimuth (true anomaly), (pj the jet polar 
angle (jet inclination=angle between the z-axis and the jet-axis) 
and / the orbital inclination (angle between the z-axis and the 
line-of-sight). The true anomaly is +90 deg at conjunctions (or- 
bital phases and 0.5). For the black hole case, Dubus et al. 
(2010) derived 39 deg and 319 deg for (pj and @j, respectively. 
Since the observer has a permanent 270 deg true anomaly (@obs) 
and 30 deg inclination ((p„bs), the jet pointing is close to the 
line-of-sight. Hence, Cygnus X-3 is more likely a "microblazar" 
rather than a microquasar. Since the true jet anomaly is some- 
what larger than that of the observer, the line-of-sight trails the 
jet during the orbit. 

Besides the jet inclination and azimuth, an additional param- 
eter in the modeling is the ratio H/d, where d is the binary separa- 
tion and H the height along the jet where high energy electrons 
are released and Compton scattering takes place. Dubus et al. 
(2010) obtained a best value of 2.66 for this ratio. As an exam- 
ple, for total masses 12 Mq and 28 Mq the binary separation d is 
3.2 Rq (2.2 xlO'' cm) and 4.2 R© (4.2 xlO" cm), respectively 
(the separation depends on total mass as M,o,'^-'). 

Numerical models for microquasar jets show a bow shock 
surrounding the jet (Perucho & Bosch-Ramon l20()8T l. We used 
two simplified geometries to represent the stationary shocked re- 
gion surrounding the jet: i) cone shells with thickness 0. Id and 
opening angles between 20-50 deg, and ii) cylinders with thick- 
ness O.ld and radius between 0.2d-0.3d and starting at distances 
between 0-1. Od from the compact star. Further, we assume that 
the shocked region immediately transforms the wind material 
into clumps of equal density with space density N,j,,„,//Nf./,„„y, ad- 
vected radially from the WR-star by velocities 1000-3000 km/s. 
N,v,„£/ and Ndump are the wind and clump densities, respectively. 
Hence, the wind mass is conserved in the process. Note that this 



Y 




Y 



Fig. 1. The system geometry (from Dubus et al. l2010l l. The jet- 
angles used are 0^ = 319 deg and (pj - 39 deg. The line of sight 
lies in a plane parallel to the yz-plane with direction {- &gi,s ) 
= 270 deg and / {-(pobs) = 30 deg. 



is just our assumption to explore the possibility of clumps with- 
out demonstrating that they really exist. In numerical computa- 
tions we used the xyz-grid divided into 70x70x70 pixels with 
pixel size 0.18d, giving 12. 6d for the box size. 

The cone geometry with opening angles less than 20 deg 
or larger than 50 deg did not produce any orbital modulation 
(Section 3), nor thinner tubes in the tube geometry. The shell 
thickness has a minimal impact since we are not interested in 
the absolute values of mass columns, but only relative values 
along the orbit. 

The geometries are shown schematically in Fig.|2]as projec- 
tions on the orbital plane. The WR-companion is in the center 
while the compact star is moving in a circular orbit around it 
(dotted circle wit four phases marked). The observation direc- 
tion is from below (see also Fig. |4]). Note that the projections 
appear to end before reaching the grid boundaries (0,70) due to 
small jet inclination (39 deg). 



3. The "clumpy trail" and the two absorbers 

We assume that the bow shock triggers the formation of dense 
clumps or makes the existing clumps in the WR wind denser. 
The formation of clumps migh physically similar to star for- 
mation when a shock front impinges upon interstellar clouds 
(van Breugel et al. 120041 1. Furthermore, we assume that these 
clumps participate in the wind outflow with velocity between 
1000-3000 km/s. Hence, once formed the clumps flow out with 
the wind. We stress that we just want to explore the possibility 
of clumps without demonstrating that they really exist. 

We assume that the clump density (like the wind density) 
decreases with distance r from the WR-star as r"'' with y-l. The 
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Fig. 2. Projection boundaries of bow shocks on the orbital xy- 
plane (pixel numbers marked) for the two geometries used at 
orbital phases 0.35 (the cone geometry) and 0.85 (the tube ge- 
ometry). The WR-companion is in the center while the compact 
star is moving in a circular orbit around it (dotted circle). The 
observation direction is from below (see also Fig.|4]l. 



Fig. 4. Clump density along the crossing of the line-of-sight and 
clumpy trail (the strip), integrated over all phases and projected 
on the orbital xy-plane (for the cone in Fig. |2]i. The WR star is 
in the center, while the compact star is moving in a circular orbit 
around it (four orbital phases are shown). The projection of the 
line-of-sight at phase 0.65 is indicated by the vertical line. 



0,00 

+++++++ 

+ + 
0.25 WR 

+++++++++ 
0.50 




20 T 




10 20 30 40 50 60 70 



Fig. 3. Projections of the shocked regions in Fig. 2 (assumed 
to be transformed into clumps) and advected by the wind 
(2000 km/s) during one hour (0.2 x orbital period). The WR- 
companion is in the center while the compact star is moving in a 
circular orbit around it (dotted circle). The observation direction 
is from below (see also Fig.|4]i. 



binary separation d = 3.4 Rq used in the computations obviously 
scales with the wind velocity. 

An example of the initial jet bow shock from many compu- 
tations is shown in Fig.|2]for a wind velocity 2000 km/s, a cone 
with opening angle 20 deg (at orbital phase 0.35), and a tube 
with radius 0.3d (at orbital phase 0.85). We further assume that 
the bow shocks transform wind material immediately into dense 
clumps, to be advected by the wind. The projections of the ad- 
vected clumps on the orbital plane are shown one hour later in 
Fig. [3] When coupling similar images over all phases and times 
we get the outflowing formation we call the "clumpy trail". Thus 
the trail is the whole volume of clumps generated by the wind 
advection. 



The crossing of the line-of-sight (from the compact star) and 
the clumpy trail was then computed numerically. For a fixed or- 
bital phase this intersection is a segment ('strip') of the line-of- 
sight along which the clump density is weighted by r^^. For the 
cone in Fig. |2] these strips integrated over all phases are shown 
in Fig. |4] as projected on the orbital plane (note the small in- 
clination of 30 deg used and twice smaller pixel size used in 
this plot). The crossing region distance r from the compact star 
is between 0.3-3d (1-10 R0) at the ionization parameter range 
log(^)=3.75-2.75 erg cm/sec for density 10'^ cm^^ {^=L/{Ni^)). 
In this plot we use a smaller pixel size than in Fig.|2](0.09d) to 
better reveal the details. 

Mass columns (assumed proportional to optical depths) 
along the line-of-sight for three cone models (20 deg, 1000 
km/s), (20 deg, 2000 km/s), (40 deg, 2000 km/s) and a tube 
with radius 0.3d are shown in Fig. [5] scaled with their maxi- 
mum values along the orbit (dotted lines). The average phase at 
maxima is around 0.35 while the FWHM is approximately 0.25. 
The solid line is from the individual cone model in Figs. 2 and 3 
(opening angle 20 deg, wind velocity 2000 km/s). Amongst the 
dozen computed models, the latter produced the best fitting re- 
sults and will be used as the clumpy trail absorber T(clump) (see 
Ch4.) 

In the same plot (Fig. |5]l we include the wind optical depth 
T(wind) (= mass column along the line-of-sight, the solid line 
with maximum around phase 0). This will be used as the first 
absorber (wind absorber T(wind)) in Section 4. Using (mass 
column)^ for the optical depth would give a better fit for the 
continuum light curve around minimum, and it may reflect that 
-y in reality is larger than 2 (accelerating wind) or else there is an 
ionization effect in the wind. Using a pure mass column yielded 
practically the same continuum light curve but with a partial 
eclipse profile that was slightly too broad. 

Around 0j = 39 deg the results are not very sensitive on 
this jet inclination angle, while changing 0^ between 290 deg- 
340 deg moves the maximum of T(clump) between 0.25-0.40 in 
phase. 
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PHASE 

Fig. 5. The two absorbers T(wind) and T(clump) used in the 
modeling in Section 4 (solid lines) where T(clump) corresponds 
in cone geometry to jet opening angle 20 deg and wind veloc- 
ity 2000 km/s. The dashed line has opening angle 20 deg and 
smaller 1000 km/s wind velocity. The two almost coincident dot- 
ted lines correspond to a wide cone (40 deg, 2000 km/s) and a 
tube (see the text in Section 3). 



The key factor in this concept is whether the clumps formed 
are radiative or not. In our case, the wind (including clumps) is 
photoionised and very radiative. Hence, the concept should work 
(from shocks to dense clumps). 

3.1. Cooling times of shocked clumps 

Bearing in mind that we have assumed, for the purpose of this 
paper, that a bow shock triggers clump formation, we neverthe- 
less estimate cooling times for clumps, both in equilibrium with 
the radiation field as well as heated to 10^ K. We speculate that 
hot clumps could exist at short distances from their origin and 
this is discussed in Section 5. 

To estimate cooling times we use the photoionisation model 
of Vilhu et al. ( 12009b . where the compact star luminosity is L ^ - 
2.46x10^^ erg/s. The clump densities probably vary, as do their 
sizes, but as an example we use a density of lO'-' cm"-'. Based on 
our XSTAR-computations with inclination 30 deg, the wind base 
density should be around lO'^ cm""' to guarantee a continuum 
optical depth around unity (at 2 keV). This is what the partial 
eclipse of the continuum light curve requires (see Section 4). A 
somewhat higher clump density is then a good guess. 

The radiation spectrum of a clump in equilibrium with den- 
sity 10'-' cm"^ at a distance of 3 R0 from the compact star 
(log(^)=3.75) is shown in Fig.|6] as computed with the XSTAR 
code (Kallman'2006). In this model, the total outward luminosity 
Lo„, per clump particle equals 4.2x10"'' erg/sec, while its ther- 
mal content per particle in equilibrium with the radiation field 
(temperature approximately 10^ K) is 1.38x10"'' erg. 

Hence, the cooling time is short (3.3 msec) but this is coun- 
terbalanced by heating from the radiation field (heating=cooling 
in equilibrium). For a clump with N = 10'^ cm"^ the cooling 
time is 20 msec (with similar equilibrium temperature 10^ K). 
Cooling times for hot 10^ K clumps (collision dominated) are al- 
most 1000 times longer: 3 sec and 17 sec for densities 10'-' cm""^ 
and 10'^ cm"', respectively. This is due to a heat content that is 
100 times larger and outward luminosities 10 times smaller of 
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Fig. 6. Radiation spectrum of a clump (erg/s/keV per particle) 
with density N - 10'"' cm"' and in equilibrium temperature 
(around 10^ K) at a distance of 3 R0 from the compact star (with 
L, = 2.46xl03** erg/s). 

the 10^ K clumps. Cooling times of 10^ K clumps can be several 
minutes. 

Hot 10^ K clumps are not in equilibrium with the surround- 
ing radiation field and try to expand (and cool) after the shock 
has passed. A crude estimate of the dynamical time scale can 
be found as follows. The force F across the clump surface area 
S is PxS, where P is the gas pressure inside a clump with mass 
m. From the clump expansion acceleration a=F/m, one can com- 
pute that the dynamical expansion time is 3.5 sec and 35 sec for 
clump sizes lO" cm and 10** cm, respectively (independent of 
clump density). 

Hence, hot 10^ K clumps can survive at most a few minutes 
after their formation, depending on their density and size. Their 
possible role is discussed in Section 5. 

4. Modeling of the X-ray continuum and emission 
line light curves 

Here we compare the observed light curves with models using 
two absorbers: i) the WR wind, and ii) an absorber lying in the 
clumpy trail (see Fig.|4|i. Both optical depths are plotted in Fig.|5] 
and explained in Section 3. 

The observed gamma-ray (ferm//LAT) light curve is taken 
from Abdo et al. (2009) using the pdf measuring perimeter tool 
for graphics (Adobe) and modeled with Dubus et al. (2010) in- 
verse Compton scattering formulas and parameters (/ = 30 deg, 
(pj=39 deg, 0y=319 deg, H/d = 2.66). The same parameter val- 
ues were used for the clumpy trail computation. Note that Dubus 
et al. (1201 Oi l used phase units where the X-ray minimum occurs 
at phase 0.25. In this paper we use the X-ray phases where the 
X-ray minimun occurs at phase (the WR star is between the 
observer and the compact star). For a jet with and Qj = 270 
deg, the jet axis and line-of-sight would coincide. 

The RXTE/ ASM light curve integrated over fifteen years, and 
scaled inside a specific orbit by the daily mean, was used as a 
template for the continuum light curve. The ephemeris of Singh 
et al. (120021 ) was used to compute the orbits. The data were lim- 
ited to moderately high states (daily means > 15 counts/sec), 
to be more consistent with the Cygnus X-3 spectral state during 
the Chandra HETG observations used. The emission line light 
curves (Si XIV Lya 6.185 A, FeXXVI Lyor 1.780 A (H-type) 
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Table 1. Fitting parameters and significances for the ASM- 
continuum (daily means > 15 cps) and Chandra emission line 
light curves. DOF is 27 and 7 for the continuum and lines, re- 
spectively. MEAN - mean of the lines. 



TYPE 


A 


B 




F-STAT 


SIGNIF 


ASM 


1.82±0.06 


0.20±0.03 


0.57 


1.00 


0.99 


SiXIV 


1.97±0.48 


1.01±0.28 


0.62 


1.18 


0.81 


FeXXV 


0.36±0.65 


0.89±0.40 


0.30 


1.25 


0.75 


FeXXVI 


1.47±0.97 


1.51±0.58 


0.38 


1.39 


0.63 


MEAN 


1.29±0.70 


1.11±0.41 


0.30 


1.23 


0.76 



Table 2. Fitting parameters and significances for the pure wind 
absorber (r(clump)=0 or B=0). 



TYPE 


A 


B 




F-STAT 


SIGNIF 


ASM 


1.50±0.04 


0+0 


2.04 


1.03 


0.93 


SiXIV 


0.55±0.34 


0±0 


2.58 


9.1 


0.003 


FeXXV 


0.0±0.0 


0±0 


1.76 


inf 


0.0 


FeXXVI 


0.0±0.0 


0±0 


1.47 


inf 


0.0 


MEAN 


0.0±0.0 


0±0 


1.40 


inf 


0.0 



''"^ FEKMI-LAT 




Fig. 7. Observed (crosses) and model (solid) light curves for 
gamma-ray (top, FermijhPil), X-ray continuum (middle, ASM 
daily mean > 15 cps) and mean of SiXIV, FeXXV and FeXXVI 
emission lines (bottom, Chandra HETG). The dotted lines in the 
two lowest panels show the contributing optical depths of the 
two absorbers (wind and clumpy trail), scaled with the maximum 
of continuum wind absorption (the same scale on the left as for 
the light curves). 

5. Discussion 

The continuum light curve is represented by the ASM aver- 
age light curve over fifteen years limited to daily means > 15 
counts/sec. So far there is no clear evidence that the form of the 
light curve changes during this time nor depends on spectral state 
- low (hard) or high (soft) - in a significant way. However, we 
limited the ASM data to relatively high states to be more consis- 
tent with the emission line observations. 

The emission line light curves were observed during a high 
state and with a few days' exposure by Chandra HETG (Vilhu 
et al. 120091 PI M. McCollough). It was assumed that they are 
representative also for the whole set of ASM observations. This 
is just an assumption and has no observational verification due 
to the limited amount of spectral data. The same applies for the 
Ferm//LAT gamma-ray observations since they are means over 
one year of observations during flaring periods (Dubus et al. 
20T0l) . 

We assumed that the eventual jet precession has a long time 
scale (over tens of years). Otherwise the observations could not 
be compared with the same jet parameters. This appears to be the 
case, since most imaging radio observations between 1991 Jan 
15 - 2001 Sept 15 locate the jet position angle in the North- 
South direction (0 deg or 180 deg; Marti et al. 2000; Mart i 
et al. I200n Miller- Jon es et a l. ^Mi Schalinski et al. [T998T l. 
However, Tudose et al. ( 12007b using 2006 May 01 observations 



and F eXXV 1.859 A (He-type) ) were taken from Vilhu et al. 
( 1200 9) and observed during a high state by Chandra/HETG (PI 
M. McCollough). The number of phase bins used was 30 and 10 
for the continuum and lines, respectively. 

The modeling of the light curves requires implementation of 
both the wind and clumpy trail absorbers explained in previous 
sections. We assume that the continuum soft X-ray emission is 
centered at the compact star (disc). Most line emission (in partic- 
ular SiXIV) comes from a broad region in the wind, except the 
high excitation FeXXVI Lya emission that probably originates 
from the compact star disc (see Vilhu et al. 2009) . Line absorp- 
tion in the continuum (at wavelengths below the line) should 
then be the common factor for all the lines. 

Let T(wind) and T(clump) be the optical depths along the 
line-of-sight (from the compact star) in the wind and clumpy 
trail, respectively, and scaled with their maximum values along 
the orbit (as shown in Fig. |5]). As baseline absorbers we adopt 
those delineated by the heavy solid lines in Fig.[5]and described 
in Section 3. 

The model flux is defined with the help of these two ab- 
sorbers as follows: 

F - exp[-AT(wind) - BT(clump)] . 

A and B are constants derived from the fitting to observations 
and given in Table 1 (using the baseline absorbers). The IDL pro- 
cedure mpcurvefit.prc^was used in the fitting. The F-statistic and 
the significance (between 0.0-1.0) of the fit was computed with 
the IDL-procedure/v' - test.pro where 1.0 represents the highest 
significance. Table 1 shows that the fits are significant for both 
the continuum and the lines. If the clumpy trail absorber is not 
used at all (see Table 2) the fits are not acceptable, in particular 
for the lines. 

Fig. |7] shows the observations overplotted with model fits. 
The mean light curve for the lines is shown, but note that Table 
1 gives fitting parameters for all the lines. The contributing opti- 
cal depths, AT(wind) and BT(clump), are also shown in the two 
lower plots by dotted lines (scaled with the continuum wind col- 
umn Ar(wind)). The scale on the left is the same as for the light 
curves. The gamma-ray {Fermi LAT) light curve is shown for 
comparison in the uppermost panel. 

It can be seen that most continuum modulation comes from 
the wind while the small asymmetry around phase 0.3-0.4 is 
caused by the clumpy trail (B much smaller than A, see Table 
1). On the contrary, emission line absorption is enhanced at the 
clumpy trail (larger B). 



' written by Craig Markwardt 
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Table 3. Integrated optical depths along the 'strip' (crossing of 
line-of-sight with the clumpy trail) at phase 0.35 (see Fig. 4 and 
text in Section 5). 



FEATURE 


ORIG WIND 


EQUIL CLUMPS 


10' K CLUMPS 


ASM 


0.04 


0.13 


0.04 


SiXIV 


0.02 


0.96 


0.14 


FeXXV 


0.12 


0.14 


1.12 


FeXXVI 


0.13 


0.05 


0.42 



locate knots (which might be interpreted as segments of jets or 
counter jets) at position angles 55 deg - 82 deg. Furthermore, 
Mioduszewski et al. (2001 ) locate the position angle during the 
1997 Feb 06 observations at 145 deg. We also assumed that the 
clumpy trail (in the circumbinary region) is more or less the same 
during all observations (including big flares or micro-flaring). 

The physical nature of the two absorbers requires more work. 
At soft X-rays below 10 keV photoelectric absorption is impor- 
tant and depends both on the distance from the X-ray source and 
wind density. It also depends on the element and line transition 
in question. Pure electron scattering (depending on the number 
of particles along the line of sight) contributes in a small way. 
To illustrate the situation we computed three different scenarios 
with the XSTAR code (using the X-ray source model explained 
in Section 3.1) along the line-of-sight from the X-ray source at 
phase 0.35 (at the strip-segment 0.3d-3d from the source, see 
Fig. 4): 

1 . Original wind - the wind was homogeneous with density de- 
creasing as r"^ with distance r from the WR-star (lO'^ cm"^ 
at the WR-surface 

2. Equilibrium clumps - the wind consisted of equilibrium 
clumps with density 10'^ cm""* (in equilibrium with the radi- 
ation field) 

3 . 1 0^ K clumps - the clumps were forced to high constant tem- 
perature 10' K 

In all cases the mass columns along the line of sight strip were 
the same (wind mass conserved). 

Continuum and line optical depths were computed along the 
strip using the XSTAR-code. Table 3 gives the integrated depths 
while Fig. 8 gives the optical depths per 10'*'cm versus distance 
from th X-ray source (in units of the binary separation). In the 
plot the clumps are presented by a mixture (50/50) of equilib- 
rium and 10' K clumps (solid lines). This mixture gives a mod- 
erate fit to the B-parameter values of Table 1, provided that the 
wind base density is twice larger 2 xlO'^ cm"^. The pure wind 
case is shown by dashed lines. The dotted lines give the contin- 
uum for clumps (upper line) and wind (lower line) cases. 

Table 3 and Fig. 8 show that the line absorption can indeed be 
enhaced at the clumpy trail. The iron lines appear to require hot 
collision dominated clumps (hotter than the equilibrium ones), 
in particular at short distances from the X-ray source. Replacing 
the clump density lO'-' cm"-' by 10'^ cm"-' or 10'"^ cm"-' does not 
change much this picture. The short cooling times (see Section 
3.1) seem to require some modification to our assumptions: e.g. 
the jet axis is closer to the line of sight, the jet cone is broader, 
the bow shock propagates with the wind or clump formation is 
delayed. 

It is noteworthy that the presence of the clumpy trail absorber 
occurs around phase 0.3 (see Fig. 3) which is also where the 9 
mHz QPOs were found by Koljonen et al. 120111 1. Whether these 




CROSSING OF LOS AND CLUMPY TRAIL AT PHASE 0.35 



DISTANCE EROM X-STAR/SEP 

Fig. 8. Optical depths vs distance from the X-ray source along 
the crossing of line-of-sight with the clumpy trail at phase 0.35 
(see Fig.4 and text in Section 5). 

QPOs arise from some sort of flickering when the line of sight 
passes the moving clumps remains to be solved. 

In the present paper a formal exponential law with two ab- 
sorbers was used to model the light curves (Section 4). However, 
it is significant to note that the location of the second absorber 
coincides with the clumpy trail region. 



6. Conclusions 

We have demonstrated a possible method to study obsei-vable 
effects of jets in microquasars when the jet direction and line-of- 
sight are close to each other (as in Cygnus X-3). In particular, 
we showed that the jet in Cygnus X-3 can produce a "clumpy 
trail" which crosses the line-of-sight during the binary orbit pe- 
riodically. 

Using jet parameters derived by Dubus et al. (2010) from 
Fermi/LAI observations, we computed the location of the 
clumpy trail (see Fig. 3). Model light curves were then con- 
structed using the two absorbers in Fig. 5 weighted with con- 
stants A and B : i) the WR wind, and ii) a clumpy trail. These 
light curves were compared with the observed ones. Good agree- 
ments were achieved for the soft X-ray continuum (RXTE/ ASM) 
and emission lines (Chandra HETG) (see Fig.|7]and Table 1). 

We found that the location of the clumpy trail computed from 
jet parameters matches well with that of the second absorber re- 
quired (clumpy trail absorber) which may not be just a coin- 
cidence. Although more work is required to clarify the physi- 
cal nature of this clumpy absorber, we suggest that the clumpy 
trail consists of a mixture of equilibrium and hot (shock heated) 
clumps. 

We note that the results are based on an assumption of con- 
stancy of the jet position angle over the past fifteen years (long 
precession time) which needs to be confirmed. 
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